Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 21 April 2009 (MN I*TeX style file v2.2) 



Matter temperature during cosmological recombination 
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ABSTRACT 

The temperature of the atomic matter in the Universe is held to that of the Cosmic 
Background radiation until decoupling at z ~ 100. After this it cools faster than the 
radiation (oc (1 + z)^ rather than (1 + z)) and would have fallen to about 20 mK 
today if astrophysical feedback processes had not heated up the intcrglactic medium. 
We show how the derivative of the Compton coupling equation helps numerically to 
follow the decoupling process. 
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At early times atoms are coupled to Cosmic Microwave 
Background (CMB) photons through Compton scattering. 
In an expanding Universe matter would 'like' to cool as 
Tm oc (1 + , i.e. faster than the radiation, which varies 
with redshift z as Tr oc + However, Compton coupling 
prevents the matter cooling this rapidly until the partial 
ionization of the atoms has fallen enough that the Compton 
heating timescale becomes long compared with the Hub- 
ble time. So although cosmological recombination is often 
referred to as 'decoupling', the huge photon bath prevents 
matter decoupling from the radiation until much later. The 
CMB 'last scattering' epoch is at 2; ~ 1100, while matter does 
not really start to cool adiabatically until z ~ 300. 

Using the WMAP Markov Chains to account for the 
variation in the cosmological model, parameters we find that 
Tm would be (0.0215 ± 0.0002) Kelvin, if there had been no 
additional sources of heat. This means that the asymptotic 
behaviour is the same as if the matter had instantaneously 
departed from the radiation at 1 + z = 2.725/0.0215 ~ 127. 
Of course the growth of non-linear structure at z < 20 and 
subsequent feedback of gravitational and nuclear energy 
leads to intergalactic medium temperatures in today's Uni- 
verse which are much higher than the CMB temperature. 

Although we found the uncertainty on today's Tm by 
considering the variation among currently acceptable cosmo- 
logical model parameters, probably a bigger uncertainty lies 

in the actual physics of recombination at low redshift. Some 
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simple algebra shows that Tm oc a;^ j , where x^j is the free 
electron fraction (normalized to hydrogen by Xe = ne/nu) 
which 'freezes out' at low redshift. The additional uncer- 
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tainty in Tm due to Xe, f is then expected to be of the order 
2-4% (see for example IChluba. Rubino-Martin fc SunvaevI 
|2003). 

The explicit equation governing the kinetic temperature 
of the matter (here meaning electrons plus ions plus atoms, 
wi th dark matter being uncoupled) is given by equation (66) 
in ISeager. Sasselov fc ScottI (|2000l ) . Ignoring the negligible 
atomic cooling processes (Bremsstrahlung, collisions, etc.) 
we have: 



(1 + 2)r^ = 



(Tm - Tr) 
Htc 



-f 2Tm, 



(1) 



where H{z) is the Hubble parameter evaluated at epoch 2 
and 

3moC 1 -f /hc + Xc 



tc 



(2) 



SaTanT^ Xc 

Here the prime denotes differentiation with respect to red- 
shift, (tt is the Thomson cross-section, qr is the radiation 
constant {— 87r^fc''/15c^/i^) and /hc is the fractional abun- 
dance of helium by number (assumed ionized here for sim- 
plicity, but correctly dealt with in the full recombination 
codes). 

Clearly Tm is slightly below Tr at early epochs, with the 
difference kept at just the right value for Compton heating 
to make the matter track the radiation. The small imbal- 
ance is also important because it produces a 'Compton drag' 
force on the matter particles. An estimate of this tempera- 
ture difference appears to have been first been mentioned by 
ICamowl l|l949l ). He states (equation 20) without proof that 

Tr - Tm f(years) 



Tm 



IQi 



(3) 



Furt her discus s ion o f this temperature difference is 
given bv lWevmaniil (|l966l h whose result is 
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Tn 



60— (1 + z) 
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(4) 



where we have converted to our notation. This has the ap- 
proximately right redshift dependence and ionization depen- 
dence, ahhough differences in assumptions about the cos- 
mological model make it difficult to compare the coefficient. 
Still, this result is essentially correct. 

The difference between matter and ra diation tempera - 
tures was also included in the textbook of lPeeblesI l|l97ll) . 
He writes a version of equation ([l} and states 'Because the 
coefficient in the last term is so very large we get a good ap- 
proximation to the solution by setting = 0' (converting 
to our notation), hence finding that 



Tn 



-H 



(5) 



This is a go od order of magn itude estimate, agreeing (in 
essence) with I WevmannI l)l966l l. but it is worth pointing out 
that TIji = would be a much better approximation than 
^1 = 0. 

Let us write Tm = Tr — e at early times, with e having 
the dimensions of temperature and fixing Tr oc [1 + z) at 
all times. Then the solution to equation ([T} is simply 



Htc 



(6) 



In the limit a;c — > 1 (and ignoring helium) this is half of the 
expression in Peebles (1971). We note that the same r esult is 
obtain ed in a rather different way in the Appendix of lHiratal 
ll2008ll . 

With this approximation in hand we can write down an 
expression for the evolution of the matter temperature by 
differentiating equation Q: 



Tr 



1 + /h 



1 + /hc + Xc 



+ 



1 + ^ 



H 



■(7) 



This expression is useful for improving the numerical ac- 
curacy of the solution to the matter temperature. The last 
two terms (in square brackets) are of similar magnitude and 
track each other at high redshift, hence can be combined. 
The second term depends on the derivative of the ionization 
fraction, and so contributes differently as a function of red- 
shift. Together the derivative in equation (0 can be used to 
evolve the matter temperature to quite high accuracy until 
the departure of Tm from Tr stops being small. 

In solving the coupled recombination equations one does 
not really need to follow Tm exp licitly at early times. In 
the commonly used code recfast llSea.ger. Sasselov fc ScottI 
1 19991 . 12OO0I : IWong. Moss fc ScottI [2OO8I ) Tm is set to Tn un- 
til Htc reaches some predefined value, and the full equa- 
tion ^ is switched on afterwards (typically at z 850). This 
leads to a 'glitch' in the solution (pointed out in lFendt et al.l 
|2008| ). We found that this glitch is easily removed by follow- 
ing equation ([7]) instead of just Tm = Tm/(1 + z) before the 
switch, and then solving the full equation ([l]) afterwards. 
The results are shown in Fig. 1. The change in ionization 
fraction, roughly 0.2% ai z = 850, leads to a ~ 0.2% correc- 
tion to the CMB power spectrum Ces when using recfast 
in a Boltzmann code. It may seem that unnecessary calcu- 
lations are being carried out by explicitly integrating the 
matter temperature at early times, but in fact the integra- 
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Figure 1. Percentage change in the matter temperature (solid 
line) and the ionization fraction (dashed line) when using the 
approximate derivative before switching to the full derivative at 
2 ~ 850. 



tor is already so fast that there is negligible effect on the 
speed at which recfast runs. 
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